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The short gravitational wave signal from the merger of compact binaries encodes a surprising 
amount of information about the strong-field dynamics of merger into frequencies accessible to 
ground-based interferometers. In this paper we describe a previously-unknown "precession" of the 
peak emission direction with time, both before and after the merger, about the total angular mo- 
mentum direction. We demonstrate the gravitational wave polarization encodes the orientation of 
this direction to the line of sight. We argue the effects of polarization can be estimated nonparamet- 
rically, directly from the gravitational wave signal as seen along one line of sight, as a slowly-varying 
feature on top of a rapidly-varying carrier. After merger, our results can be interpreted as a co- 
herent excitation of quasinormal modes of different angular orders, a superposition which naturally 
"processes" and modulates the line-of-sight amplitude. Recent analytic calculations have arrived 
at a similar geometric interpretation. We suspect the line-of-sight polarization content will be a 
convenient observable with which to define new high-precision tests of general relativity using gravi- 
tational waves. Additionally, as the nonlinear merger process seeds the initial coherent perturbation, 
we speculate the amplitude of this effect provides a new probe of the strong-field dynamics during 
merger. To demonstrate the ubiquity of the effects we describe, we summarize the post-merger 
evolution of 104 generic precessing binary mergers. Finally, we provide estimates for the detectable 
impacts of precession on the waveforms from high-mass sources. These expressions may identify 
new precessing binary parameters whose waveforms are dissimilar from the existing sample. 



I. INTRODUCTION 

Coalescing comparable-mass black hole binaries are 
among the most likely and useful sources of gravitational 
waves for existing and planned gravitational wave de- 
tectors like LIGO [JJ, Virgo, 0, the Einstein telescope 
[3], and proposed space-based detectors. For sources in 
a suitable mass range, the signal these detectors receive 
contains significant features from the late-stage, strong- 
field dynamics of the black hole merger. Only full nu- 
merical simulations of Einstein's equations can provide 
first-principles models for this epoch, including all dy- 
namics and emission [JJ [5]. Given the large computa- 
tional cost per simulation, relatively few well-determined 
models have been produced. 1 

For sufficiently high mass mergers, ground-based grav- 
itational wave detectors are sensitive only to the short 
portion of the waveform provided by numerical simula- 
tions. Historical comparisons of the leading-order emis- 
sion from comparable-mass nonprecessing binaries has 
suggested these waveforms look similar [§]■ Exist- 
ing upper limits on gravitational waves from > 50Mq 
merging binaries imply astrophysically plausible signal 
amplitudes for the first few detections must be low, at 
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1 Most simulations have thoroughly explored nonspinning and 
spin-aligned binaries |5j. Some simulations with more generic 
spin have been performed (e.g., to study gravitational wave re- 
coil kicks) [6HH]. 



best near the detection limit of the advanced LIGO and 
Virgo detectors |10| . Short, similar, low amplitude sig- 
nals can carry a limited amount of information. Nonethe- 
less, strong-field dynamics produce significant evolution 
in the emission beampattern and polarization content 
during merger. If in band, these features should be ac- 
cessible to gravitational wave detectors. 

In this paper we provide the first phenomenological 
description of the "kinematics of merger" . Rather than 
discuss properties of horizon or of the global spacetime 
[TT1 IT!?] , we emphasize phenomenology of the asymptotic 
radiation. In particular, we identify precession-induced 
changes in the polarization content as an easily-accessible 
phenomenological measurable with which to distinguish 
between simulations. While we adopt a simple model to 
characterize this diagnostic, our results can be translated 
to fully-developed parameter estimation and model selec- 
tion strategies, which fit data to a model using physically- 
motivated [T5HT5] and phenomenological [T5H2T] signal 
parameters. For clarity we use three examples to illus- 
trate the precession we have in mind and the imprint it 
makes on emitted signals. To demonstrate the features 
we describe occur frequently and to connect those fea- 
tures to quasinormal mode frequencies, we also describe 
the post-merger dynamics of a much larger array of sim- 
ulations. In Section|ll]we describe the gravitational wave 
signals, preferred orientations, and description of polar- 
ization we use in our work. We argue our practical ap- 
proach to gravitational wave polarization, decomposing 
the line of sight signal into left- and right-handed com- 
ponents, is both observationally accessible and intuitive. 
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Our study uses detailed numerical relativity simulations 
of the late inspiral and merger of two binary black holes, 
in general both with spin. By example, we demonstrate 
that "precession" dynamics, suitably defined, continues 
after merger. Again using selected examples, we illus- 
trate that each line of sight encodes the relative orien- 
tation of the (time dependent, precessing) preferred ori- 
entation. By comparing these to a simple expression, we 
argue this behavior occurs in general, for all lines of sight 
and all precessing binaries. In Section |Hl| we demonstrate 
that all these features change significantly, depending on 
the spin magnitudes and orientations adopted. Based on 
these investigations, we anticipate the polarization con- 
tent can trace features of the strong-field merger event. 
Then, in Section |IV[ we demonstrate that this polar- 
ization information is experimentally accessible: similar 
waveforms differing primarily through their polarization 
content (i.e., time history of left- versus right-handed) 
can be observationally distinguished. We demonstrate 
observations can weakly constrain orientation evolution 
nonparametrically , independent of other information 
naturally encoded in the gravitational wave signal. 



A. Context 

Though the gravitational- wave signal from merging bi- 
naries nominally encodes all information about the space- 
time, previous investigations of the (2, 2) mode suggest 
universality of the merger signal [5J [§] , providing the in- 
formal impression that previously hypothesized compli- 
cated nonlinearities in general relativity do not signifi- 
cantly complicate the short signal from the merger epoch. 
In other words, simulations suggest the merger signal pro- 
vides little unambiguous information about strong-field 
general relativity. However, the unique contribution of 
merger to the waveforms is invariably mixed with the 
preceding (inspiral) and subsequent (ringdown) signal, 
making nonlinear behavior difficult to identify. In princi- 
ple the mass loss, angular momentum change, and recoil 
kick can all be tabulated, fit, and constrain the nonlinear 
response at merger. Nevertheless, in practice the contri- 
bution from merger cannot be cleanly distinguished from 
the contributions from inspiral and ringdown. While 
symmetry demands a general form |22j . symmetry does 
not strongly forbid inspiral or merger epochs from con- 
tributing to these integrated quantities. The contribution 
from merger can be partially disentangled in the time 
domain, plotting fluxes like dE/dt, or possibly the time- 
angular domain. Unfortunately, investigations looking 
for strong mode-mode couplings indicative of nonlinear- 
ity have recovered the familiar perturbative results (231 — 
|2"5] , Previous investigations suggest that the asymptotic 
radiation encodes preferred directions with different sym- 
metry properties, not tied to the orientations associated 
with standard conserved constants. These orientations 
are well resolved through merger and seem to evolve non- 
trivially both before and after the merger [57] . This 



preferred orientation may therefore provide a new signa- 
ture of strong-field dynamics. 

After this study was completed, we became aware of 
an analytic study providing a simple, geometric inter- 
pretation of (short-wavelength) quasinormal modes' fre- 
quencies as precession [28]. Their calculation extended 
previous estimates derived in the slow-rotation limit to 
arbitrary spin. We find a similar result at the lowest mul- 
tipolar order, even in the strong field: all modes precess 
at a nearly constant rate during and after merger. 



II. SIMULATIONS AND DIAGNOSTICS 
A. Simulations I: Overview 

The simulations examined in this paper will be de- 
scribed at greater length in a companion publication. 
Physically, these binary configurations correspond to an 
unequal-mass generalization of the "S series" [pj, where 
one spin is in the orbital plane ai = ax and the other 
spin rotates through the (x, z) plane with comparable di- 
mensionless spin a = |S/Af 2 |. Initial data was evolved 
with Maya, which was used in previous binary black 
hole (BBH) studies |>H3 [HH33] ■ The grid structure for 
each run consisted of 10 levels of refinement provided by 
CARPET [33], a mesh refinement package for CACTUS [55] , 
Sixth-order spatial finite differencing was used with the 
BSSN equations implemented with Kranc |36) . Initial 
binary properties and simulation details including initial 
separation, spin configuration, and finest resolution are 
detailed in Table [I] In the text, we present results for 
the fiducial simulation at r = 90M unless otherwise in- 
dicated. The choice of extraction radius has minimal 
impact on the results shown here. 

To consistently illustrate features of these simula- 
tions, we will emphasize one fiducial binary, denoted 
Sq(4,0.6,90,9) with m 1 /m 2 = 4, a% = 0.6a; = — a 2 , placed 
in orbit in the x, y plane wth a coordinate separation 
d = 9M; see Table [i] for details. 2 Figure [l] shows the 
corresponding inspiral-merger trajectories of the orbital 
angular momentum and spin directions, L(t) and Si(t), 
where S±(t) corresponds to the larger black hole. Ini- 
tially, the binary's orbital angular momentum is larger 
but comparable to Si = aim 2 ; both are much larger 
than the companion's spin. Early in the inspiral, L(t) and 
Si (t) precess around the total angular momentum direc- 
tion [37j . During this process, the three angular momenta 
(J, L and Si) are nearly coplanar, with the angle between 
L and Si roughly constant, as expected from the post- 
Newtonian limit. By contrast, as the orbit transitions 
from inspiral to plunge, the (coordinate) orbital angular 



2 Our fiducial initial conditions are comparable to those adopted 
in a companion paper on orientation-dependent emission |27l . 
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FIG. 1: Evolution of angular momenta and V with 
time: For one simulation (Tq(4,0.6,90), with mi/m2 = 4, 
ai — 0.6x, a2 = 0.65 starting at d — 10M), a plot of the 
preferred direction V (thick black; see Eq. superimposed 
on the coordinate orbital (black) and spin (green) angular 
momenta directions L, Si as a function of time. The paths 
indicate how all three evolve with time; the arrows are evalu- 
ated at a specific instant shortly before merger. A blue arrow 
also indicates the final (~ initial) total angular momentum 
J. At very early times, the directions V and L nearly agree 
|38| ; at merger, however, L converges to J while V continues 
to precess. Compare to figures in O'Shaughnessy et al. |27] . 



momentum converges towards the total angular momen- 
tum direction. The spin direction, by contrast, changes 
far less dramatically during the plunge and merger. 



B. Preferred orientations 

Conserved quantities of the spacetime and radiation 
escaping to infinity encode several preferred orienta- 
tions. One preferred, manifestly physical, but nearly- 
time-indenpendent orientation is J, the nearly conserved 
orientation of the total angular momentum. Another ori- 
entation extremizes the quadrupolar radiation, such that 
the (2, 2) mode of the Weyl scalar is largest (3^1 BD]; we 
denote this orientation as Q(t). Finally, a third orienta- 
tion V(t) is the principal axis of [2"rJl |2"?1 13"5] 



(£(„£&)) 



Jdfl\i^\ 2 

Ezmm' falm'i'ilm ( lm ' |A>A) | lm) 



(1) 



/ dfl\^ 4 \ 2 

where C are rotation group generators; where "04 is the 
Weyl scalar; and where in the second line we expand 



V^4 = £/ OT V>4i m (*) y jL (9,<f>) and perform the angular 
integral; see Ochsner and O'Shaughnessy [38] for explicit 
formulae. 

These three quantities do not agree. For example, early 
in the inspiral J is manifestly distinct from the two pro- 
cessing, dynamic orientations V, Q, which both nearly 
correspond to L [38, 39]. The two dynamic orientations 
also generally differ, albeit less significantly. 3 As a con- 
crete example, if the Weyl scalar consists only of the I — 2 
subspace with the time-independent value 
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then on the one hand our preferred direction has V ~ 
n(0.05, 0) but the direction that maximizes the (2, 2) 
mode lies along n(0.25, 0), where n(6,<p) is a cartesian 
unit vector with polar coordinates 8, <j). For simplicity, 
in this work we only present results for a single dynamic 
orientation V. 

Figure [l] shows the trajectories of L, V and J. While 
L and V nearly agree prior to merger, they differ signifi- 
cantly at and after merger. In the time domain, the direc- 
tion V roughly corresponds to the direction of peak | ^4 1 • 
More critically, when averaged over mass and frequency, 
this direction corresponds very closely to the direction of 
largest signal amplitude 
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where we adopt the isotropic two-detector network def- 
inition of p 2 used in O'Shaughnessy et al. [57]. Mov- 
ing forward in time, both V and L precess around the 
total angular momentum J [37] . When the two black 
holes merge, the orbital angular momentum has con- 
verged to the total angular momentum (L — > J). This 
direction no longer evolves past merger. By contrast, the 
observationally-relevant direction V continues to evolve 
through and beyond merger. 

As will be demonstrated using the fiducial simulation, 
different preferred directions identify different physics. 
In particular, V is directly connected to observations, 
and indicates interesting new dynamics that could be im- 
printed on the gravitational wave signal. 

Here, we identified this imprint at each time using 
all gravitational wave emission directions. Experiments 
have access to only one line of sight. Below, we describe 
a concrete method to identify the imprint of V using only 
a single emission direction, via the polarization evolution 
versus time. 



3 While the differences between the two approaches to preferred 
orientations are often small in absolute scale, these small an- 
gular differences imply significant disagreement about the rela- 
tive magnitude of strongly suppressed modes, particularly in the 
merger phase. 
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C. Line of sight diagnostics: amplitude, phase, 
polarization 

Real gravitational wave detectors only have access to 
a single line of sight. Despite this limitation, they still 
have access to the rich and distinctive modulations that 
the orbit and precession imparts on the line-of-sight time- 
evolving gravitational wave signal, h = h+ + ih x . To in- 
terpret these modulations, we change basis from the stan- 
dard linear polarizations (h + ,h x ) to two circular polar- 
izations (hfi, Hl with h = h R + hi 1 ). Strictly, we split the 
Weyl scalar ip4 rather than h(t), to minimize errors from 
poorly-constrained early-time and late-time effects. 4 We 
will see below that precession distinctively modulates the 
amplitude and phase of both circular polarizations. 

The left and right-handed circularly polarized parts 
of ipi(t) are naturally isolated in the frequency domain 
[27j : positive frequency components are right-handed, 
and negative frequency components are left-handed. In 
the time domain, this procedure corresponds to convolu- 
tion of h with a specific (acausal) kernel: 
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where Kr^l(t) are suitable inverse fourier transforms of 
the unit step function. For each polarization, we can de- 
fine an amplitude \hfj\ and phase arg(Zi^). In general, 
the amplitude evolves on a precession (or, if none, radia- 
tion reaction) timescale; the phase increases secularly, on 
a radiation reaction timescale, with precession-induced 
modulations on shorter timescales |37j . 

This procedure for separating left- and right-handed 
signals can be perfomed with real gravitational wave 
data. For each line of sight, experiments with comparable 
sensitivity to two linear polarizations can reorganize their 
data analysis procedure to be sensitive to only these left- 
or right-handed polarizations. Additionally, this projec- 
tion process has physically expected properties, applied 



4 We anticipate the gravitational wave content of the (2, 2) mode 
to be almost exclusively right-handed. However, the linear trans- 
formation from h — > fiR L is highly nonlocal, depending in prin- 
ciple on arbitrarily early and late times. The gravitational wave 
strain at the start (inspiral) or end (nonlinear memory) can be 
significant. By contrast, the Weyl scalar's diminished early-time 
amplitude and lack of memory terms leads to a well-behaved 
projection. 



to nonprecessing binaries. The angular modes ip&i m de- 
scribe emission preferentially above (m > 0) or below 
(to < 0) the orbital plane, where to is the mode order 

(i.e., each term in i/>4 = 2Zi m ^HmXim ^)' s mode de- 
composition is proportional to some "04/m expim4>). One 
can empirically verify that the angular modes from non- 
precessing binaries are nearly chiral: 



fpAlm(f) - for to/ < 



(5) 



In the stationary-phase limit, this relationship corre- 
sponds to a desirable and usually satisfied requirement on 
the angular frequency versus time: 5 9targ-04; m is mono- 
tonically increasing for m > and decreasing for to < 0. 
For nonprecessing binaries, to a good approximation our 
polarization projection corresponds to eliminating modes 
of either to < (projecting to R) or to > (projecting 
to L). 

To decouple polarization content from secular trends 
in amplitude and phase seen in nonprecessing binaries, 
we once again change variables from ip^ R L (2 complex 
or 4 real quantities) to a real typical amplitude A, a real 
typical phase $, and a complex polarization amplitude z 
(one complex and two real quantities): 
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By construction, the polarization amplitude z is nearly 
constant for nonprecessing binaries whose gravitational 
wave emission is dominated by equal in magnitude and 
conjugate in phase (l,m) — (2,2) and (2,-2) modes. In 
this limit, the Weyl scalar asymptotically takes the form 

rMn,t) = M [e i*<t)y£-a> ( ft) + e -«(t)^ (fl)] (g) 



To recover the correct limits along the z axis (i.e., the ±z 
direction) , the basis coefficients of the two spin- weighted 

(—2) 

harmonics Kj-ta are necessarily right and left circularly 
polarized. As a result, for a signal dominated by equal 
and conjugate (2, 2) and (2, —2) modes, the polarization 
amplitude is a time-independent and purely geometrical 
expression: 
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where 9 = cos -1 L-h and \&£ characterize the orientation 
of the orbital angular momentum along the line of sight 
and in the plane of the sky, respectively; see, e.g. [4"T] . 



Due to projection effects of spheroidal harmonics onto spheri- 
cal harmonics, this property need not hold for all higher-order 
modes. We will not address angular mode mixing in this paper. 
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Reversing the sense of the approximation above, as- 
suming instantaneous and symmetric (2, ±2) emission 
along a proposed slowly- varying orientation 0(t) , we 
estimate the polarization content Bo(t) along a line of 
sight h by 



zo = 



[6-{x n + iy n )] A {l-d-n) 2 
(l-(O-n) 2 ) 2 (1 + 6 -h) 2 

[6 ■ (gn + ^n)] 4 
(1 + 6-n) 4 
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where the two vectors x n , y n define a frame of reference 
in the plane of the sky, perpendicular to n. As this sim- 
ple approximation suggests, for a precessing binary the 
polarization amplitude should fluctuate, reflecting the 
relative importance and phasing of left- and right-handed 
emission in the instantaneous signal. 

Closer to the orbital plane, the balance of polariza- 
tions is nearly equal (\z\ ~ 1). This naturally finely- 
tuned region involves near-perfect cancellation of (some 
of) the leading-order emission, 6 allowing higher-order an- 
gular dependence to contribute significantly to the polar- 
ization content. Unlike the leading-order factors, these 
higher-order multipoles sometimes treat the component 
masses asymmetrically. As a result, the waveform and 
polarization content near the orbital plane has additional 
modulations at the orbital period. For example, for a 
nonprecessing binary with both conjugate and symmet- 
ric emission in all modes, the Weyl scalar has the form 



(12) 



As instantaneous quadrupole emission almost always 
dominates, the polarization amplitude z can be approx- 
imated by an oscillating but nearly-constant correction 
to the leading order term [Eq. ( 10 )] : 
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In short, for orientations where higher harmonics con- 
tribute significantly, the polarization amplitude z^ 
should oscillate, with the peak-to-trough amplitude in 
one-to-one relation to the magnitude of these harmonics. 

To summarize, any gravitational wave signal can be de- 
composed into two polarizations ip^ R ^. Using the com- 
plex amplitude z = ipi^/i^iR to characterize the relative 
proportion and phasing of each signal, we find that po- 
larization content encodes information about the merging 
binary on at least two scales. On the one hand, on long 



precession timescales, the polarization amplitude tells us 
how the preferred emission direction evolves; see, e.g., 
Eq. (11). On these scales, the complex amplitude can 



6 Near the oribital plane of a nonprecessing binary, gravitational 
wave emission is nearly linearly polarized. 



very significantly. On the other hand, on the orbital and 
eventually merger timescale, rapid small oscillations in 
the polarization amplitude directly measure the relative 
proportion of (2, ±1) modes. These modes are produced 
naturally in asymmetric or eccentric binaries along cer- 
tain lines of sight. Finally, in our analysis we have for 
simplicity assumed symmetric emission in the (2, ±2) 
modes. However, when expressed in a corotating frame 
[26] . many of our simulations emit asymmetrically, with 
one mode preferentially larger than another during the 
merger and ringdown phase. Without foreknowledge of 
this bias, the simple estimates used above would, if in- 
verted, recover the incorrect inclination of the preferred 
emission direction relative to the line of sight. In other 
words, while the complex polarization amplitude z man- 
ifestly encodes information about the preferred orienta- 
tion's evolution relative to the line of sight and about 
the relative significance of higher modes, we cannot un- 
ambiguously interpret this information unassisted by a 
large catalog of candidate waveforms. For the purposes 
of this paper, we will treat z(t) itself as a phenomenolog- 
ical observable. 



D. Line of sight diagnostics: Three examples 

To illustrate the power of this decomposition, Figures 
[2] [3] and [4] apply it to a nonprecessing and precessing 
q — 4 binary, respectively. 

Figure [2] demonstrates our polarization decomposition 
with a nonprecessing q = 4 binary, extracting the Weyl 
scalar along a "generic" direction [(9,4>) — (60°, 205°)]. 
For a nonprecessing binary, the gravitational wave signal 
is dominated by quadrupole emission from the (2, ±2) 
modes. To an excellent approximation, the two polariza- 
tions ifriRL (red and blue curves in this figure) are pro- 
portional to one or the other of these modes, respectively, 
with the proportionality constants set geometrically by 
spin- weighted harmonics [Eq. ([9])] and hence the incli- 
nation. However, this line of sight is close enough to 
the orbital plane that the (2, ±1) mode can contribute 
significantly. As demonstrated with the dotted line in 
the bottom panel of Figure [2j these modes beat against 
the leading-order quadrupole, causing the complex po- 
larization z to oscillate by roughly ten percent during 
the inspiral at roughly the orbital period. Higher-order 
angular modes I > 2 also beat against the leading- 
order quadrupole. Including these terms in the the sum 
Sjm ^imXim produces richer time dependence along 
each line of sight, in both polarizations (solid curves in 
the top panel of Figure [2]). Nonetheless, the ratio of the 
two polarizations (z) still behaves like the / = 2 result. 

Precessing binaries, by contrast, generally exhibit dra- 
matic changes in polarization. Figure [3] shows features 
of the precessing mx/m-i = q — 4 binary started with 
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FIG. 2: Polarization imprints on signal 1: Nonprecess- 
ing: For a nonspinning q = 4 binary, two measures of the 
polarization content along (9,<f>) = (60°, 205°). Results are 
shown including all angular modes ip4 lrn with I < 4 (solid) and 
just retaining terms with I — 2 (dotted). Top panel: The line 
of sight right and left-handed amplitudes \ip4 R L \ (red R, blue 
L). For this and almost all other lines of sight, a single helicity 
dominates for all time. For comparison, the solid black line 
shows \ip4\ = \lp4, R +ip4 L \. Bottom panel: Kez^, = Reip4* L /tp4 R 
vesus time. Though one helicity dominates, the polarization 
content oscillates significantly at the orbital period. These 
oscillations are proportional to the "higher-mode" content 
(2,±1) and therefore measure the mass ratio or residual ec- 
centricity. 



a\ = 0.6a; = —02 at a coordinate separation d = 9M, as 
used in [27]. As seen in the top panel, the balance be- 
tween the two polarizations changes significantly: though 
R-handed emission (red) usually dominates, for a short 
epoch L-handcd emission is stronger in this direction 
(blue). Critically, we can identify when this transition 
occurs by comparing our line of sight with the preferred 
orientation V(t) extracted from (£( a £ b )) (251 [57] • Fol- 
lowing the discussion above, we expect and our calcula- 
tions confirm (shaded region) that the two polarizations 
have equal amplitudes (i.e., \z\ = 1) when V ■ h = 0, 
with R-handed emission dominating when V ■ h > and 
L-handcd emission dominating when V ■ n < 0. More 
precisely, we can estimate both the magnitude and phase 
of z surprisingly reliably by combining Eq. ( 11 ) with the 



FIG. 3: Polarization imprints on signal 2: Process- 
ing: For a precessing binary, the polarization carries a strong 
imprint from the relative orientation of the preferred emis- 
sion direction relative to the line of sight. For the precessing 
q = 4 binary described in the text and |27| . a decomposition 
of the gravitational wave signal along a "generic" orientation 
[(8, (j>) = (60°, 205°)] is shown. Results are shown using just 
/ = 2 (solid) and all modes I < 4 (dotted), reversing the 
convention of Figure [2] Top panel: The line of sight right 
and left-handed amplitudes \ip4 R L \ (red R, blue L) using all 
modes (dotted) and just I = 2 (solid). The gray shaded box 
shows the interval where V ■ ri < 0. Bottom panel: Com- 
parison of \z(t)\ extracted along this line of sight (black solid 
[/ = 2] and dotted curves [I < 4]) with the leading-order esti- 



mate z^y provided by Eq. (Ill and the preferred orientation 



orientation V(t). For example, the bottom panel of Fig- 
ure [3] shows the simulated polarization amplitude \z\ is 



V selected by (£( a £(,)) (blue) 



extremely close to the estimate \zy\- Small fluctuations 
around the leading-order prediction at the orbital period 
should be due to residual eccentricity. 

Despite emphasizing I = 2 modes in the discussion 
above, the inclusion of higher-order angular modes I > 2 
into ip4 produces nearly no change to \z\. In this ex- 
ample, our naive expression relating z to a preferred di- 
rection continues to apply, even when these harmonics 
produce dramatic fluctuations in the line-of-sight ampli- 
tude. Judging from this example, similar calculations for 
the phase of z, and their repeated success for all lines 
of sight and simulations considered, we believe z is pri- 
marily determined by the orientation of fi relative to the 
time-dependent preferred orientation V and can usually 



be well- approximated by Eq. (11) 
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As described in subsequent sections, our calculations 
suggest the preferred orientation V evolves significantly 
and rapidly during merger. Equivalently, the polarization 
content - the distribution of orientations dominated by 
left versus right handed emission - changes significantly 
at the merger event. As seen in Figure |3J immediately 
before and after the merger event, the two points cor- 
responding to predominantly left- or right-handed emis- 
sion change noticably. This interval corresponds to the 
merger phase itself. As described in the next section, we 
suspect this rapid, global change in polarization content 
may reflect features of the strong- field merger event itself. 



E. Waveforms along other fixed directions 

For reference, in Figure[5]we show the polarization con- 
tent for two other preferred orientations: the initial (~ fi- 
nal) total angular momentum direction J (top panel) and 
the preferred orientation V evaluated at the time of peak 
emission (bottom panel). In the first case, one polariza- 
tion is vastly larger than the other at early times; during 
the merger, however, both polarizations become signif- 
icant. Similar results are found when extracting along 
z, the initial orbital angular momentum. In the second 
case, both polarizations are comparatively large early on. 
During the merger epoch, however, only one polarization 
dominates. Generally speaking, when adopting a fixed 
frame one can choose to simplify some narrow epoch of 
the waveform by reducing the other polarization. For 
any time, frequency, or mass range, a generalization of 
(£( a £f,)) can be constructed to determine what orien- 
tation would be suitable. However, in general no one 
orientation works for all time. 

Finally, we emphasize that we have been able to accu- 
rately estimate the polarization content using the time- 
dependent preferred orientation V . For our simulations, 
this orientation differs substantially from L at all times. 
Based on this performance, we anticipate that corotating- 
frame waveforms along V will be substantially simpler 
than any analog extracted along L. 



III. SIMULATIONS II: TRENDS AND 
VARIATIONS 



From the diagnostics above, we anticipate simulations 
are best and naturally characterized by (a) the modal 
waveforms ip4[ m m a corotating frame and (b) the evolu- 
tion of our preferred orientation with time. In addition, 
to simplify the translation between time and frequency 
domain, we will also use (c) the overall orientation- 
averged signal power p. In this section we briefly report 
on salient ways these three features change with spin and 
mass ratio. 



A. Polarization bias 

Generic precessing binaries exhibit a polarization bias: 
at any given instant, the binary is radiating more of one 
handedness than another. During the inspiral, the bal- 
ance between L and R oscillates. At merger, the balance 
fixes, preferring one handedness, with the choice depend- 
ing on the spin-orbit configuration just prior to merger. 
As described in Healy et al (in prep), this asymmetry 
produces large kicks, with a significant component per- 
pendicular to the orbital plane; cf. [32] . 

For the purposes of this paper, this asymmetry com- 
plicates our interpretation of the preferred direction. As 
with nonprecessing binaries, we use the ratio of left to 
right-handed power to estimate an (instantaneous) incli- 
nation. The expressions that we invert for O ■ n [e.g., Eq. 



(Ill] assume an equal amount of left and right handed 



power. More refined estimates that relate the line of sight 
to the left/right ratio are required when a strong bias to- 
wards one or the other handedess occurs. 

To assess whether more complicated expressions would 
be required in parameter estimation, 7 we define polarized 
analogs of the orientation-averaged signal amplitude p 
g3], starting from p 2 [Eq. H)]: 



- P \h) = 

! 2 



dn df \Mf,n)\ 2 

4tt S h (2tt/)4 
df \Mf,n)\ 2 



dn 

-in 

4^ 2 Jo S h (27T/)' 
dft c 

4-7T ' 



df \MfM 

S h (2tt/)4 



(14) 
(15) 
(16) 



Any given binary has directions where one polarization 
dominates (e.g, along L for a nonprecessing binary). 
Summing over all orientations, however, nonprecessing 
binaries emit symmetrically, with matching amounts of 
R and L handed emission along mirror-symmetry-related 
lines of sight. With equal amounts of signal power, non- 
precessing binaries must have Pl = Pr = P/V%- By 
contrast, precessing binaries have no symmetry that en- 
forces symmetric emission; particularly in narrow epochs 
or frequency intervals selected by an outside observer or 
gravitational wave detector, one handedness (R or L) can 
dominate. As Figure [6] demonstrates with a specific ex- 
ample and Figure [7] with an ensemble, generally a single 
polarization does dominate during merger. The domi- 
nant polarization depends sensitively on the spins, par- 
ticularly when Si , S2 are nearly antiparallel and in the 
orbital plane. 



7 We illustrate polarization asymmetry using a detection-weighted 
diagnostic to demonstrate the bias is detectable. Alternatively 
and in a detector-agnostic way, the polarization asymmetry also 
shows up clearly in (J dfl\^f>ji(n)\ 2 )/(J dfl\tp£ J (n)\' 2 ), where the 
numerator and denominator are evaluated at each time or each 
frequency. 
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FIG. 4: Evolution of amplitude, polarization at merger: Three snapshots of \ip4\ (top panel) and \z\ (bottom panels) for 
the Sq(4,0.6,90,9) simulation, bracketing the time of peak amplitude and demonstrating that the polarization content changes 
significantly during the merger event. For aesthetic reasons, we only show the contributions from all / = 2 modes. 



Our calculations suggest that ~ 30% changes in the 
overall amplitude (p or ^4) are not uncommon during 
merger. The relative amplitude of left to right handed 
radition (pr/ pl or z), however, changes by orders of mag- 
nitude due to small changes in inclination; see Eq. (11) 
and Figure [3j As a result, the typical polarization bias 
introduces a fairly small systematic error into any pro- 
cedure to reconstruct the evolution of V [e.g., into the 
inverse of Eq. (11) for the inclination as a function of 
time] . 



modes are chiral in a way that is relevant to data analy- 
sis. For example, we can decompose the contribution of 
each corotating mode to ficticious 8 "orientation-averaged 
signal-to-noise ratios" pr,l associated with each handed- 
ness: 



Pirn ,R,corot 



Phn,L,corot 



df i«r(/)i 2 



Si, 



47T* 

1 ( 



df WJ£ 



(/)l 2 



(17) 
(18) 



B. Corotating waveforms, chirality 

The corotating-frame waveforms will be described in 
detail in a subsequent publication. For the purposes of 
this paper, we will employ only two key features. First 
and foremost, like their nonprecessing analogs, the coro- 
tating frames are chiral: modes with m > have fre- 
quency content only for / > and vice-versa. We 
can quantify how precisely we are confident that the 



Applying these expressions, we find Pi m ,L/pim — for 
m > : modes that should have positive helicity (m > 0) 
have little negative-frequency power. As a corollary, the 
corotating-frame (2, 2) and (2, —2) modes are nearly or- 
thogonal - as, more generally, are the subspaces spanned 



Ficticious since they are associated with a corotating frame. 
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FIG. 5: Polarizaton content along fixed special direc- 
tions: The right- and left-handed Weyl scalar amplitude 
\ip4 RL \ (red, blue) versus time, for the gravitational wave 
signal extracted along several preferred orientations. For 
comparison, the gray lines show the simulation-frame (2, ±2) 
(solid gray) and (2,±1) modes (dotted gray). For aesthetic 
reasons, we only show the contributions from the I = 2 modes. 
Top panel: Initial total angular momentum. Bottom panel: 
V at the time of peak emission. From the ratio of left to right 
(blue to red), the direction identified by V corresponds to 
nearly circular polarization near the epoch of peak emission. 
This direction is significantly offset from J. 



by modes with m > and m < 0. We will use the orthog- 
onality of the corotating-frame (2, 2) and (2, —2) modes 
to model how precisely we can measure the orientation 
of the corotating frame, both on average and as function 
of time. 




FIG. 6: Merger signals are polarized overall: Depending 
on the spins, binary black holes can emit substantially more 
L or R handed power. As an example, a comparison of pR 
(red) and pL (blue) to p/v2 (dotted) for a binary starting 
with ai = an(60°),a2 = az. In this expression we adopt a 
fiducial initial LIGO design noise curve. 
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C. Preferred direction precesses 

For the simulations and time intervals we have simu- 
lated, the preferred orientation V evolves as if precessing 
along a nearly-constant cone, centered along some axis 
W. Figure ^ provides a concrete example. In almost all 
cases we find W ~ Jfinal empirically: V precesses around 
the total angular momentum. For this section only, how- 
ever, we allow W to take arbitrary values. 

Given that orientation and an arbitrary frame xw, Vw 
defined perpendicular to the constant vector W, we de- 
fine the precession phase (j)w and precession frequency 



PlIPr 



FIG. 7: Large L/R asymmetries and large kicks: Scat- 
ter plot of the recoil kick velocity in the z direction (top panel) 
and overall (bottom panel) versus Pl/pr- Large kicks in the 
+5 direction correlate with predominantly L-handed emis- 
sion. The largest possible kicks, however, correspond to fairly 
moderate asymmetries between L and R during merger. In 
this figure, colors correspond to different simulation sets: S 
(blue), Sq (red), T (yellow), Tq (green). For illustrating the 
correlation between R vs L bias and kick magnitude only: kick 
data is computed using the simulation series alone at a single 
extraction radius, without correcting for the early inspiral. 
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FIG. 8: Precession of preferred orientation: Exam- 
ple: Demonstration that the preferred orientation precesses 
around W — J for q = 2, d = 10A/, ai — 0.6a;, ai = 0.6z. 
Top panel: Two-dimensional path of V (blue) and L (green) 
around a proposed value of W (point at origin), as seen in 
a plane perpendicular to W. The two dashed curves are 
circles of constant opening angle. For comparison, a blue 
point indicates the final angular momentum direction J in 
this frame. Center panel: J ■ V (blue) and J ■ L (green), 
showing both L and V precess along similar cones, with a 
nearly constant opening angle prior to merger. In this fig- 
ure, merger occurs at t ~ 0. Bottom panel: Plot of <j>w before 
and after merger, demonstrating an abrupt change in the pre- 
cession frequency at the merger event. For comparison, the 
green curve shows (j>L, the precession phase extracted from L 
around J. For this system, at late times the "precession rate" 
ilvw = dt<j>w — 1/12M is still an order of magnitude smaller 
than the "carrier frequency" set by the (2, 2) and (2, —2) 
modes. 



\(x w + iy w ) ■ V\ 
flvw = dt<pw 



(19a) 
(19b) 



and make similar definitions for L around J. During 
the inspiral, the two vectors V and L nearly coincide, 
precessing around a similar center at a similar rate. To 
quantify this similarity, Figure [8] shows the opening angle 
of each precession cone (J ■ L and J ■ V) and the "phase" 
4>v.J an d 4>lj of the vectors L and V as they precess 
around J. Prior to merger, both quantities largely agree. 

As the binary plunges and merges, the precession cones 
of both L and V cycle around their respective axes (J and 
W) at a higher rate. By comparison to L, however, the 
preferred axis V has a much smaller precession rate and 
persists in precession cycles long after merger. Further- 
more, contrary to our intuition, the opening angle of the 
precession cone for V (i.e., arccos(VF • V)) rarely decays 
significantly after merger. This precession cycle of V is 
robust: both for this simulation and all others tested, all 
extraction radii and (when available) all resolutions show 
quantitatively similar features. 

For each simulation, we estimated W and 4>w by fitting 
a fixed precession cone to the trajectory of V. Specifi- 
cally, we adopt as W the direction that mimimizes the 
rms difference between V(t) ■ W and its time average. 
This choice corresponds to the assumption that the pre- 
cession cone swept out by V should have a constant open- 
ing angle. When estimating W, we separately fit data 
prior to merger and after merger. In almost all cases we 
found W ~ J and flyj — Slvw- 

To further demonstrate the close correspondence be- 
tween these two directions, Figure [9] shows both an- 
gular frequencies: points appear at the median value 
(Qvw + ^V./)/2; the bar has height ±(Qvw + £V,/)/2. 
Except for a handful of cases, mostly associated with 
nearly-aligned spins, these two expressions agree. We 
conservatively adopt the difference (£lyw~ ^Vj)/2 as an 
estimate of our systematic error. [Other sources of sys- 
tematic error, such as extraction radius, have a smaller 
effect on the recovered post-merger precession rate.] 



D. Post-merger precession 1: Results 

In our simulations, if the binary precesses prior to 
merger, our preferred direction V continues to pre- 
cess after merger, never converging towards J. Preces- 
sion is ubiquitous (i.e., whenever L and J are initially 
misaligned); we therefore expect ubiquitous precession- 
induced modulations at and after merger when asym- 
metric binaries merge with randomly-oriented spins. As 
described in a companion publication, we have performed 
extensive resolution tests and remain confident that the 
post-merger oscillations we observe are resolved. Table [T] 
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provides an estimate of the post-merger "precession fre- 
quencies" fiyj for each of our simulations. 

Figure [9] demonstrates that the post-merger "preces- 
sion" frequencies are principally determined by the fi- 
nal black hole's spin magnitude. Because the final black 
hole can be described as a superposition of quasinormal 
modes, the precession frequencies are very nearly differ- 
ences between the (l,m) — (2,2) and (2,1) quasinormal 
modes of the remnant hole [44] : see Figure [9] 

By some change of basis, the long-lived post-merger 
oscillations must be described as some perturbation of 
the final remnant black hole: a (coherent) superposition 
of multiple quasinormal modes. Our calculations suggest 
precessing mergers coherently excite multiple quasinor- 
mal modes of different angular order. These quasinormal 
modes have similar decay timescales (1/Im(w)) but not- 
icably different characteristic frequencies (1/Re(w)). As 
a result, though in principle the quasinormal modes are 
not exactly degenerate, in practice just after merger these 
coherently-excited quasinormal modes decay at a similar 
rate but with changing relative phase, leading to a black 
hole state that appears to precess. 

To illustrate why our preferred direction precesses, we 
consider a toy model that approximates features of a su- 
perposition of outgoing quasinormal modes: 




a2\m\e 



ity(-2) 
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FIG. 9: Post-merger precession frequency: For each of 
the simulations listed in Table [T] a scatterplot of their "post- 
merger precession frequency" Af/flvw [Eq. (19 1] against the 
final black hole's spin magnitude Jf/Mf. Colors indicate four 
different classes of simulation: Tq (blue), T (red), Sq (green) 
and S (yellow). For comparison, the solid lines show the 
quasinormal mode frequency differences ^122— ^121 (blue) and 
UJ222 — W221 (yellow), where ujnim is the real part of the quasi- 
normal eigenfrequency [44] . The bars indicate differences be- 
tween two methods for estimating Q,; see Section [HI C| 



E. Post-merger precession 2: Broader context 



where cj; m is the lowest-order eigenfrequency for the 
(l,m) eigenspace, where the eigenmodes are defined rel- 
ative to the black hole's angular momentum direction 
J, and where a 2 \ m \ are three real parameters. In this 
example, we ignore the difference between the actual an- 
gular eigenf unctions of the Kerr-background wave equa- 
tions (spin-weighted spheroidal harmonics) and the corre- 
sponding eigenfunctions for flat space, Y^ n ' . In the triv- 
ial case where only a 2 ±2 are nonzero, the preferred orien- 
tation remains along the black hole's angular momentum 
axis J |45j . Allowing a 2 ±i to be nonzero but keeping 
020 = leads to a preferred orientation V that spirals in- 
ward exponentially towards J, precessing around z with 
a frequency f2y = Re(w22 — W21) an d shrinking towards z 
at a different rate Ty = Im(uj 22 — ^21)- Finally, when all 
mode amplitudes are nonzero, the preferred orienation 
can exhibit a wide range of behaviors depending on the 
relative mode ratios (121/(122 and 020/022, including ex- 
ponential decay to J; modulated precession around J; 
and pathologically complicated behavior in the presence 
of degeneracies. 

In this interpretation, the angle between V and J re- 
flects the amplitude of the (2, ±1) modes. Conversely, 
the angle between V and J is also partially geometric, 
being tied to the spin-orbit configuration and how it the 
binary evolves during merger. However, as is apparent 
from Figure [sj the late-time (2, 1) amplitude (i.e., the 
angle V ■ J) evolves during the strong-field merger event. 
We anticipate the opening angle provides new insight into 
the strong-field merger process. 



Our results suggest that the nonlinear merger event 
coherently seeds multiple quasinormal modes. After 
merger, all quasinormal modes appear to precess almost 
coherently, in a frame aligned with the final black hole's 
spin. For example, each constant-/ subspace seems to 
precess at a similar rate. We defer a detailed discussion 
of each constant-/ subspace to subsequent publications. 

Our results agree with a recently-developed geomet- 
ric interpretation of Kerr quasinormal modes. Using a 
short-wavelength limit, Yang et al. [H] recently showed 
that these kinds of quasinormal mode frequency differ- 
ences are intimately connected to precession of spheri- 
cal photon orbits. Our calculations suggest this kind of 
geometric interpretation holds even in the strong field, 
during merger. 

Because rapid polarization changes occur during the 
merger epoch, we anticipate that these polarization 
changes encode linear and potentially nonlinear features 
of strong-field gravity. At a minimum, the precession 
cycle of V persists during and after the merger, with a 
precession rate tied to the quasinormal mode frequency 
distribution of the final hole. Modulations in the line 
of sight polarization content provide a natural way to 
identify this characteristic frequency, potentially allow- 
ing new experimental tests of general relativity. Second, 
the long-term, multimodal coherence we observe suggests 
a short, nonlinear, and nontrivial process seeds a coher- 
ent multimodal perturbation of the final black hole. In 
part, of course, the perturbed black hole simply inher- 
its features from the just-prior-to-merger binary. In the 
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TABLE I: Simulations used: The first column is a key, encoding the family, mass ratio, black hole spin magnitude \Si\/M 2 — |S 2 I/ M 2 anc ^ alignment. 
The next 8 columns provide specific initial conditions: the initial separation (r s t ar t) : mass ratio g = m-i/m^, and two component spins S 2 / M 2 relative to the total 
initial mass. The next column provides P£,/pn, a measure of whether the merger event preferentially radiates R or L-handed emission. [This quantity is evaluated for 
M — 200Mp| and the initial LIGO design noise curve.] The next column provides £ly j , the post-merger "precession frequency" of V around J. estimated using C± 50M 
after peak (2, 2) emission; "unk" indicates entries where we cannot reliably determine it. The last three columns provide the simulation duration, the length of the reliable 
waveform, and the highest resolution h used. 
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geometric point of view, V ■ J 7^ prior to merger; in 
the quasinormal mode view, the (2, 1) mode just prior to 
plunge partially limits the (2, 1) quasinormal mode am- 
plitude during and after merger. However, the merger 
process partially but incompletely sheds this quantity. 
We anticipate that the amplitude of this effect, measured 
either by V ■ J or the peak (2, 1) amplitude in a frame 
aligned with J, will provide a useful probe of the strong- 
field merger process. 



F. Recovering the direction using polarization 



Detailed Fisher matrix and Markov-chain Monte Carlo 
calculations corroborate this simple estimate [T71 146H52] . 

For nonprecessing binaries dominated by I = \m\ = 2 
emission, the polarization content is equally a measure 
of the relative amplitude (and phase) of the (2, 2) ver- 
sus (2,-2) modes; see Eq. |l0| ). Because the (2,2) and 
(2,-2) modes have definite (and opposite) helicity, we 
can equally interpret inclination measurements and po- 
larization content constraints as information about the 
relative importance of left- and right-handed emission, 
averaged over the whole signal. 



For processing binaries, we have argued that the line- 
of-sight polarization, as a function of time, determines 
both the line of sight and the path of V, a preferred 
direction. We have argued this process works for all time, 
both prior to and after merger. In the theorists' paradise 
- access to a noise-free line-of-sight signal ip4(t) - we have 
already arrived at this paper's key result: by tracking 
polarization, we can partially reconstruct properties of 
the source binary, even when limited to the short merger 
epoch. 

In practice, however, gravitational wave detectors have 
limited sensitivity, with access to only a small time and 
frequency interval of the full signal. Nontheless, as we 
argue below, both the polarization bias, pl/pr, and pre- 
ferred orientation, V, can be measured by tracking the 
time-dependent polarization along a single line of sight. 



IV. POLARIZATION IS MEASURABLE 

In the previous section we described how to charac- 
terize the polarization content of a gravitational wave 
signal. In this section we demonstrate that this con- 
tent is experimentally accessible, both on average and 
to a lesser extent as a function of time. We quantify the 
amount of polarization difference that experiments can 
distinguish. In particular, we show how fairly small po- 
larization fluctuations can lead to substantial differences 
between nonprecessing and precessing waveforms. 



B. Static and dynamic polarization content 
accessible 

As demonstrated extensively for nonprecessing signals, 
gravitational wave observations can constrain the polar- 
ization content of a static signal. By implication, ob- 
servations can equally well constrain the "static part" 
of a signal with weakly-time-varying polarization con- 
tent. Generalizing, observations should also be able to 
constrain some "average polarization" from an arbitrary 
gravitational wave signal. How? Conceptually, gravi- 
tational wave detector networks can be arranged to be 
sensitive to only one (circular or linear) polarization at 
a time. By measuring the total power incident in each 
polarization, we directly constrain the "average polariza- 
tion." 

A detailed analysis of the response of real gravitational 
wave detector networks to each circular polarization is 
beyond the scope of this paper. Instead, for simplicity 
and following the philosophy outlined in |27j and near 
Eq. m, we will idealize gravitational wave networks as 
equally sensitive to both linear polarizations. By con- 
struction, the left- and right-handed components of ^4 
are orthogonal. The net signal amplitude p 2 can be ex- 
pressed as a sum of left- and right-handed components: 



4 



= 2 



df \Mf,n)\ 2 
S h (2tt/)4 

df \Mf,n)\ 2 
,S h (2tt/) 4 



(21) 
(22) 



A. Nonprecessing binaries 



As shown above [Eq. (10)], for nonprecessing binaries 



measurements of the polarization content z are equiva- 
lent to constraints on the emission inclination (9) . While 
phrased in different coordinates and applied to an en- 
tirely different mass regime where semianalytic waveform 
models exist, several extensive discussions of parameter 
estimation including inclination exist in the literature 
[PTl 25-52 . For a signal of amplitude p, the inclination 
can be measured to order 

A(cos6l) ~ 1/p 



Observations determine the integrated contributions 
from both polarizations (pr.l) independently. In this 
averaged sense, observations can unambiguously deter- 
mine the polarization content: we can simply measure 
the (source mass and detector-dependent) ratio 



z(M) = PL I Pi 



(23) 



For nonprecessing binaries, the average polarization am- 
plitude corresponds precisely to the instantaneous polar- 
ization amplitude given in Eq. (10). Motivated by the 



accuracy of inclination measurements and fluctuations in 
each (independent) polarization, we anticipate this quan- 
tity can be measured to order Alnz ~ l/min( j o^, p^). 
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For sufficiently strong signals, the Weyl scalar can be 
reconstructed in the time domain, allowing us to esti- 
mate ip4R l an d therefore the polarization content at each 
time. Nonparametric signal reconstruction algorithms 
like coherent waveburst are already employed in at- 
tempts to identify and reconstruct signals from unmod- 
eled sources [531 HI] ■ At the signal amplitudes expected 
from the first merger detections, however, fully non- 
parametric models generally have too much freedom to 
tightly constrain the polarization evolution seen in merg- 
ing black hole binaries. Nonetheless, tighter constraints 
should be achievable, to the extent that precession- 
related timescales remain long compared to other scales 
in the signal. 

In fact, several post-Newtonian parameter estima- 
tion studies have demonstrated how well these longer 
timescales can be constrained and differentiated, using 
both ground- and space-based interferometers [49j l5T| 
152] . Implicitly, the initial conditions of an inspiralling, 
precessing binary encode its precession trajectory, albeit 
in suboptimal coordinates. 9 To the extent that post- 
Newtonian study of binary parameter estimation con- 
firm both spins and directions of merging binaries can 
be measured, they also imply the (time-dependent) pre- 
cession trajectory can be distinguished. Furthermore, 
post-Newtonian simulations suggest that even marginal 
precession (i.e., even less than one precession cycle) has a 
dramatic impact: if precession is included, then precess- 
ing signals can appear dramatically different than their 
similar nonprecessing counterparts. 



C. Polarization-induced mismatch: Loss of 
amplitude relative to nonprecessing 

To demonstrate the critical impact even weak but time- 
dependent precession has on a merger signal, we remove 
it. Specifically, we start with our standard precessing bi- 
nary, examine all possible emission directions, and deter- 
mine how well a comparable nonprecessing analog could 
fit them. To reduce ambiguity, for our "nonprecessing" 
analog we use the corotating-frame (2, ±2) modes from 
that same simulation. 10 The procedure for constructing 
these waveforms is described in a companion publication 



9 The traditional coordinate system for binary parameters adopts 
the initial conditions and spin vectors, rather than the geometry 
of the precessing binary as it passes through the sensitive band. 
Recent calculations suggest coordinates adapted to the center of 
the band and which phenomenologically encode precession will 
more transparently represent the available information; see, for 
example, Brown et al. |41| . 
10 Alternatively, one can estimate how much precession matters by 
comparing a precessing NR signal with a nonprecessing template 
family, including a full search over all possible component masses 
and emission directions. This challenge will be addressed in a 
subsequent publication, using existing and recently developed 
numerical and theoretical signal models. 



|55) . As noted above, these two corotating modes are 
orthogonal, with support only for / > or / < respec- 
tively. Specifically, a nonprecessing analog of our binary 
would be well-approximated by any linear combination 
of the (2, 2) and (2, -2) modes. 

As in our prior work |26|, I27j , for simplicity we com- 
pare two complex waveforms using an inner product that 
accounts for both polarizations simultaneously: 



(2tt/) 4 ^ 



(24) 



In this expression we are not maximizing over time or 
phase. Using this inner product, we can define two nor- 
malized L and R-handed basis states from the corotating- 
frame modes: 



\ R ) = ^422(0/v / (' ( /'422>422) 
\L) = ^42-a(*)/\/ (^42-2 > ^42-2) 



(25) 
(26) 



Using these orthonormal basis states, we can construct 
the "nonprecessing signal" ip es t that best resembles 
ip(t,h): 

\i> est ) = (R,1i)\R) + (L,il>)\L)=Vil> (27) 

where V is a projection operator to the subspace spanned 
by R, L. From the difference between this state and the 
original state, we can determine how much information is 
lost by a nonprecessing approximation for each candidate 
line of sight: 



2 

Prern 



^,(l-V^)=p 2 -(\(R,^\ 2 + \(L^)\m 



This expression can be further simplified by expanding 
if) in left- and right-handed functions of time along that 
line of sight: 



= PL V'L^+Pfl i>R) = 4>}p 



(29) 



where for clarity we adopt hats to denote normalized sig- 
nals (e.g., tf\ — \ip) I p). Substituting this expression 

then using orthogonality of L and R handed signals gives 
an expression for the SNR lost in terms of the intrinsic L 
and R handed amplitudes pb.,l and the overlaps (i?, ipft) 
and {L,Tpi,) between the two (normalized) line of sight 
R and L basis signals and our two R- and and L-handed 
basis functions: 



Prem 



= 1 - 



1 



Pi 



pi 



[pl(\(Ri>R)f 



pI\(lM)\ 2 ) 

(30a) 



In the above discussion we pessimistically perform a 
point-to-point comparison: we do not allow the R, L 
waveforms to vary in phase or time, neither relative to the 
signal ip nor relative to one another. In practice, however, 
we want to mimic the response of a template manifold, 
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including all physical signals that include our synthetic 
nonprecessing binary. Allowing for similar events at dif- 
ferent times, a better diagnostic for lost signal power due 
to precession is 



that any nonprecessing signal cannot mimic sufficiently 
strong geometric precession-induced phase and ampli- 
tude modulations. 



rem,mm 



mm t p r 



(30b) 



where the minimum is over all candidate times t for the 
merger event specified by \R) ,\L). 

This "missing power" diagnostic behaves qualitatively 
differently for nonprecessing and precessing binaries. For 
simplicity, consider a scenario dominated by just the 
leading-order (/,m) = (2, ±2) emission dominates, so 
higher harmonics can be omitted to a first approxima- 
tion, such as M = IOOMq and the initial LIGO noise 
curve. For a nonprecessing binary and for masses domi- 
nated by quadrupole emission, the two overlaps \ (R, iPr)\ 
and \(L,ip L )\ are nearly unity. Hence, almost no sig- 
nal power is lost. By contrast, precession introduces 
distinctive phase modulations that cannot be produced 
by nonprecessing sources [40l [41] . Since precession oc- 
curs during merger, we expect loss of signal amplitude 
in direct proportion to the phase and amplitude modula- 
tions that precession introduced. In the language of the 
above diagnostic, we expect that for precessing binaries 
the two overlaps \ (R, ipn)\ and \ (L,i/jl)\ are substantially 
less than unity. 

As an example, Figure |l0| shows contour plots of these 
two overlaps (top panel) and of the total signal ampli- 
tude lost (bottom two panels) for a fiducial precessing 
binary. In this figure, we adopt a low reference mass 
(AI = 100M Q ) to insure both that I = 2 emission dom- 
inates and occurs along a relatively stationary instanta- 
neous direction; see, e.g., Figure 2 in [27]- Even in this 
limit, the bottom left panel of Figure 10 shows that most 
orientations are poorly fit with an optimal nonprecess- 
ing approximation generated from the corotating frame. 
Further, these poorly-fit orientations lie nearly perpen- 
dicular to some instantaneous V direction. Finally, in 
each of these poorly-fit orientations, the relative phase 
of L and R changes significantly. As seen in the bot- 
tom right panel, while the phase of i\)^ R L individually 
resemble a nonprecessing waveform, together their rel- 
ative phase evolves in a way we expect that no single 
nonprecessing system could reproduce. Similar behavior 
occurs for all the precessing binaries we have explored. 

Though we adopt a particularly well-chosen reference 
signal, this calculation only concretely demonstrates a 
this particular precessing binary has many lines of sight 
where its signal cannot be fit by this particular pair of 
nonprecessing basis signals. Real nonprecessing searches 
compare each individual signal with a template bank of 
many candidate nonprecessing signals. We will address 
this "fitting factor" problem in a subsequent publication 
that compares our precessing signals to a family of non- 
precessing sources. However, motivated by other stud- 
ies of precession- induced modulations [41] , we anticipate 



D. Estimating the mismatch 

Motivated by low-mass results [41] . our calculations 
suggest a familiar conclusion: when viewed near its in- 
stantaneous orbital plane, a precessing binary can pro- 
duce complicated emission, delicately balancing left and 
right handed power. No nonprecessing signal model can 
reproduce it. Of course, our signal model also neglects 
higher harmonics. Some of the mismatch seen in Fig- 
ure [? ] reflects our inability to recover these other 
modes. For the mass region considered, however, the 
modes not included produce only a few percent of the 
(corotating- frame) signal power; see Table 1 in [27]. The 
lack of higher harmonics cannot explain the extremely 
low matches seen here. 

Low-mass studies suggest precession introduces coher- 
ent phase and amplitude modulations, superimposed on 
top of secular (corotating- frame) evolution 0T]. A suf- 
ficiently generic nonprecessing model can fit any secu- 
lar phase, but not the modulations. To show this, we 
will rewrite the line-of-sight left- and right-handed signal 
ipi R L (t) as a modulation factor times a corotating- frame 
signal. Doing so will let us calculate the overlap between 
corotating-frame modes and any line of sight as a sim- 
ple integral, whose integrand has two factors: (a) secular 
terms, reflecting the corotating modes' slow change in 
amplitude with time; and (b) a geometrical, precession 
term X, reflecting the change in amplitude and phase 
along the line of sight. These expressions prove the signal 
power "lost" in comparing a nonprecessing to precessing 
model can be uniquely associated to the modulations that 
precession introduces in the signal model - modulations 
that no nonprecessing signal can reproduce. For brevity 
and for practical reasons - almost all search strategies 
use the (2, ±2) modes alone - we emphasize a sufficiently 
low reference mass that the signal is dominated by the 
(I, m) = (2, 2) and (2, —2) modes in its corotating frame. 
As noted above, higher order harmonics / > 2 are still 
significant at these masses [Table 1 in 27 \]. Nonetheless, 
we explicitly neglect contributions from I > 2, both to 
demonstrate our analytic control over the problem and 
to connect to nonprecessing searches. 

The line-of-sight gravitational wave signal ip4(n,t) can 
be expressed in terms of corotating-frame harmonics 
i^ilm" an d a rotation operation: 



ripi(h, t) 



E D 2 mm ,(R(m 



,Z ot Yt 2 \n) (31) 



where R(t) is a suitable rotation (e.g., the forward trans- 
formation z — > V(t)) and D l mm , is a Wigner D matrix 
describing the action of SU(2) rotations on angular eigen- 
states. If the corotating-frame (2, ±2) modes dominate, 
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FIG. 10: Nonprecessing approximation omits signal power: Comparison of corotating (2, ±2) subspace with the sum of 
all Z = 2 modes along each line of sight, for the Sq(4,0.6,90,9) simulation and M = IOOMq. This binary preferentially emits 
L-handed signal into our detector's sensitive band [Pl/pr = 0.76 from Table [I]. Top left panel: Contours of normalized matches 
\(R,^Pr)\ (white) and \(L,tf>i)\ (black) for the Advanced LIGO noise curve at lOOMpc. For comparison, the colors indicate 
proximity to a single direction (h(6, (j>) ■ V* for V* = n(1.01, 2.98)), the estimated preferred emission direction at 1OOM0. The 
R basis signal is a good match when it dominates emission and vice-versa. Top right panel: Contour plot of p 2 em>m i n / p 2 , the 
fraction of signal power lost in a nonprecessing approximation [Eq. ( 30 1] . [This calculation allows independent timeshifts in L 
and R and therefore underestimates the true fraction of power lost.] In directions nearly perpendicular to the preferred direction 
V* defined above, a nonprecessing approximation fails to capture all available signal information. Bottom left panel: Fraction 
P (<) of orientations with p 2 em ,min/ p 2 l ess than a specific threshold (i.e., which lose less than a specified fraction of the signal 
amplitude with a natural nonprecessing approximation). A significant fraction of all orientations are significantly impacted 
by precession. Bottom right panel: For the line of sight (6, cj>) — (0.94,0.94), the phase of ip4 R (magenta) and ip4 L (blue) are 
shown. The two differ by a significant, time-varying phase (arg(z)). While each polarization still resembles a nonprecessing 
signal, no single nonprecessing signal model can fit both R and L polarizations simultaneously. For this reason, along this line 
of sight a large fraction of the signal power is lost when fitting with a nonprecessing approximation (top right panel). 



the line-of-sight signal can be approximated by just two 
terms: 



the rotation operator R (i.e., on V). 



r^ est (n,t) ~ ^ c 2 ° 2 rot XR + ^?X L (32) 

X R = J2 D m2(R(t))Y 2 { - 2 \n) (33) 

m 

X L EE J2 D m-2(R(t))Yt 2 \n) (34) 



The two slowly-varying but not constant functions Xr l 
depend only on geometric factors: our line of sight h and 



The variation of X^l with time is completely responsi- 

we 



ble for the structure seen in Figure 10 In Sec. IV C 



described a nonprecessing-signal search strategy, where 
each line of sight was projected into the span of two 
corotating- frame basis signals cx tp4,2±2 ■ K the two 
modulation factors Xrl were constant, these projections 
would be trivial; the search would recover the line-of- 
sight signal. To estimate the fraction of signal power lost 
due to precession, we further approximate the right- and 
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12345 12345 

FIG. 11: The line of sight estimate, Eq. |32[ is a good approximation for the full 1 = 2 waveform: Comparison of 
ip4 with the sum of all I — 2 modes along each line of sight, for the Sq(4,0.6,90,9) simulation and M = 1OOM0. As before, the 
colors indicate proximity to a single direction (n(9,4>) ■ V* for V* = n(1.01, 2.98)), the estimated preferred emission direction 
at WOMq. Left panel: Contour of the normalized match, | (tp4, eat , t[>4)\, for the Advanced LIGO noise curve at lOOMpc. Right 
panel: Contours of normalized matches \(R, 4> est )\ = 1(7?,^')! (white) and \(L,ip est )\ = | (L, ipl st )\ (black) for the Advanced 
LIGO noise curve at lOOMpc. For comparison, the thin dashed lines are the same as those in Fig. |10[ s top left panel. 



left-handed parts of ip 4 (h,t) [Eq. Ml] by 



(35) 
(36) 



which is an excellent approximation to the extent that 
Xjt i are slowly-varying. Given this ansatz, the formulae 
described in Sec. IIV CI let us calculate the fraction of 
signal-to-noise p rem / p that a nonprecessing search should 
miss when applied to this particular line of sight. For 
example, the overlap (R, i^r) can be approximated by 
(i?, tpji ), which can in turn be approximated as 

(R\X R \R) 



{R\W\R) X ' 2 
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(37) 
(38) 

(39) 



where we have simplified the integral limits using chiral- 
ity. A similar expression applies to the other chirality: 

(L\X L \L) 
121 r\V 2 



{L\\X L \>\LY 



(40) 
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f) 4 s h 



Combined these two approximations allow us to approx- 



imate Fig. 10 These approximations suggest that mis- 



match results from fluctuations in X around their median 
values. To be concrete, the total SNR lost due to apply- 
ing a nonprecessing approximation can be expanded as 



2 
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(42) 



For a nonprecessing binary with a fixed preferred direc- 
tion, the geometrical quantities X remain constant. By 
contrast, for a precessing binary, X fluctuates, with the 
fraction of lost power being proportional to rms fluctua- 
tions in X. 

Though simple, these expressions are surprisingly pow- 
erful. In particular, they allow us to estimate the critical 
threshold for precession, above which its effect is obser- 
vationally accessible, at a given mass and amplitude. As 
a simple rule of thumb, a physical process or parameter 
can impact a Bayesian posterior if the match changes by 
of order 1/ p 2 for p ~ 10 the signal amplitude. Using non- 
precessing waveforms i/'2±2 and synthetic precession tra- 
jectories V, we can efficiently explore what lines of sight 
and what precession trajectories permit an observation- 
ally significant amount of precession. We will return to 
this subject in a subsequent publication. 

As in previous studies of precession |41) . our analy- 
sis lends itself to geometric interpretation and approx- 
imation. For example, early in the inspiral, the or- 
bital frequency is both well-defined and shorter than all 
other scales. A separation of timescales argument sug- 
gests that the time-domain multiplication Xr \ R) could 
be implemented as a frequency-domain multipolication 



18 



Xnx{t{f))- In this limit, the inner products (R,^n) 
and (L,iPl) become (frequency) averages over the sig- 
nal. This approach has been used to calculate fitting fac- 
tors between nonprecessing searches and precessing black 
hole/neutron star binaries 0T]. Further investigation is 
needed to determine whether a separation-of-timescales 
approximation remains useful during the merger and 
ringdown phase. 



V. CONCLUSIONS 

Our paper makes two claims: first, that significant 
"orientation" or "polarization" changes occur after the 
merger of two binary black holes; and, second, that these 
changes could be accessible even to present-day detectors, 
if a source could be detected. 

First and foremost, in this paper we report on and pro- 
vide a simple phenomenological interpretation of rapid 
changes in the gravitational wave emission of merging, 
precessing black hole binaries at and beyond merger. 
Phenomenologically, the features we report on corre- 
spond to coherent, multimodal oscillations in the merger 
and ringdown signal that resemble "precession" of the 
remnant black hole. We confirm their existence through 
two time-domain methods: first with "polarization" , the 
ratio of left- to right-handed emission; and second using a 
preferred direction that traces the direction of strongest 
emission. We provide a simple rule to explain the ampli- 
tude and phase of these oscillations in terms of coherent, 
long-lived "precession" of the black hole binary before 
and the merger remnant after merger. Our interpreta- 
tion agrees with a recently-developed geometric interpre- 
tation for high-frequency Kerr quasinormal modes |28j . 
This correspondence merits further investigation. 

Second, using data-analysis-motivated diagnostics, we 
demonstrate each line of sight carries information about 
its orientation relative to the preferred orientation de- 
scribed above. For the simulations we have performed, 
we are able to explicitly show that a nonprecessing search 
will perform poorly, in extreme but not uncommon cases 
losing ~ 20% of signal power even at low mass (M ~ 
100M Q ). For sufficiently loud signals, we furthermore 
suggest that the polarization can be measured nonpara- 
metrically. By analogy with low mass binaries, we ex- 
pect that the most rapid phase evolution characterizes 
the binary component masses and spins, while secular 
modulations encode the orientation and spins. To sum- 
marize, the short signal from merging waveforms carries 
a surprising amount of information about merger dynam- 
ics. Our analysis suggests this information can be phe- 
nomenologically extracted, even without a complete and 
physically-parameterized model for the waveform. 

Our study has significant implications for present data 
analysis and parameter estimation strategies. First and 
foremost, our proof-of-concept study shows how data 
analysis strategies can extract easily-understood physics 
from complicated, short, multimodal gravitational wave 



signals. Second, our analysis suggests that nonprecessing 
models (phenomenological; ringdown; et cetera) for high- 
mass sources (M > 100M Q ) can omit significant features 
and introduce uncontrolled biases when applied to "typ- 
ical" merger events, as typical events can involve signif- 
icant post-merger oscillations. For example, our stud- 
ies suggest that typical moderate-spin mergers produce 
a coherent, multimodal ringdown signal, with significant 
polarization changes along a typical line of sight. As an- 
other example, by analogy to studies of inspiralling low 
mass binaries |41j , nonprecessing searches will be particu- 
larly ineffective for certain lines of sight and spin configu- 
rations. We will address the selection bias of nonspinning 
searches for precessing binary mergers in a subsequent 
publication. 

From our preliminary investigations, we believe the 
surprisingly rich and rapid behavior of merging, precess- 
ing binary black holes may encode nonlinear features of 
strong-field gravity. For example, our analysis suggests 
that rapid orientation changes occur at merger, encod- 
ing a precession frequency. We anticipate that the pre- 
to post-merger transition coherently seeds a variety of 
quasinormal modes, with surprisingly long-lived coher- 
ence in the time domain. We suggest the relationship 
between the pre- and post-merger amplitudes of these 
coherent post-merger oscillations will provide valuable in- 
formation about the strong-field merger process. All of 
the features described above can be extracted from sim- 
ulations and observations using the techniques described 
in this work. We therefore suspect that strong polar- 
ization changes are observationally accessible features of 
strong-field gravity during binary black hole merger. 
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